Coulomb Blockade in a Coupled Nanomechanical Electron Shuttle 
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We demonstrate single electron shuttling through two coupled nanomechanical pendula. The 
pendula are realized as nanopillars etched out of the semiconductor substrate. Coulomb blockade 
is found at room temperature, allowing metrological applications. By controlling the mechanical 
shuttling frequency we are able to validate the different regimes of electron shuttling. 
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One of the fundamental experiments in classical mechanics is the coupled pendula, used to demonstrate how energy 
is transferred from one resonator to another. This coupled pendula or coupled oscillator model is often evoked as an 
analogy for the Josephson effect in superconducting junctions, where similar differential equations can be applied [T]. 
fn principle, the physics of coupled oscillators can be found in mode- locked lasers [2|, analyzing weather fronts in 
climate modeling [3], and neural networks 

Nanotechnology delivers the tools to fabricate nanoscale electro-mechanical systems, pushing towards the ultimate 
limits of miniaturization |SHS]. This has already led to the realization of single nanomechanical pendula in various 
forms, e.g. suspended semiconductor cantilevers p!0Hl2] and nanopillars [TSHIS] . Such a mechanical resonator is 
commonly placed between two electrodes, so that the pendulum can exchange electrons and mechanically transfer 
them from source to drain. Hence, the term electron shuttle was coined [121 [13 • Here we present measurements of two 
coupled pendula realized as nanopillars integrated between two contacts. The strongly reduced co-tunneling in two 
serially coupled electron shuttles reveals Coulomb blockade (CB) at room temperature. The fundamental importance 
of CB in nanomechanical shuttles is two-fold: one in metrology application as outlined by Weiss and Zwerger |18| and 
the second in ultra-sensitive nanomechanical sensors. 

In classical electron turnstiles [19j, the stochastic nature of tunneling is typically suppressed by operating the 
turnstile at frequencies /, much lower than the inverse of the time constant r (or, more formally, the 1/t = 1/RC- 
frequency of the contacts). In contrast to this, we are able to operate electron shuttles in the regime of high-frequencies 
{i.e. f ^ T"^), and show that clocking of electron transport can be achieved. This effect is due to the suppression 
of co-tunneling when the shuttles are operated in series [H]. The two shuttles considered here are placed between 
two nano-scale contacts. We characterize the device's response by probing the direct current through the nanopillars. 
Mechanical motion is studied first by coupling a DC-bias and then, exciting the mechanical motion by adding a radio 
frequency (RF) signal to the source electrode, as we have demonstrated before pO] . 

The two nanopillars are defined on a silicon-on-insulator (SOI) substrate, where the top crystalline silicon is 190 nm 
thin and the insulating Si02 is about 350 nm. A 50 nm top gold layer serves as the final electrical conduction path. 
The deposited metal is also employed as a mask in a dry etch step, which mills out the SOI material around the 
pillars. We apply a CF4 plasma etch step and mill into the Si02 insulating layer, thus ensuring electron transport 
via the metallic islands. Further details on the processing are given by Kim et al. [14]. In Fig. 1(a), the final sample 
is shown in a scanning electron microscope graph. The source and drain contacts are placed in close proximity to the 
two pillars. The distance between the pillars is about 17 nm. The gating electrode is placed further away, enabling a 
shift in the electrostatic potential of the islands. The inset of Fig. 1(a) gives a broader view of the coplanar- waveguides 
into which the two nanopillars are embedded. 

The nanoscale circuit is placed in an impedance matched transmission line in order to minimize signal loss along 
the line. All measurements are performed under vacuum (< 10~^ mbar) in a probe station at room temperature. The 
station is placed in a Faraday cage and equipped with radio frequency contact probes covering the range from DC 
to 50 GHz (a bias-tee allows AC/DC superposition with high precision). The equivalent circuit diagram is given in 
the inset of Fig. 1(b): The two pillars are individually displaced by Xi and X2, which lead to tunable resistances and 
mutual capacitances (arrow boxes). The gating electrode couples capacitively to the pillars. The output current of 
the coupled electron shuttles is fed into a current amplifier. 
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FIG. 1: (a) Two coupled electron shuttles realized as nanopillars. The metallic top layer allows electron exchange with the 
source (S) and drain (D) contacts. The scale bar corresponds to a length of 200 nm. The inset shows a broader view of the 
coplanar-waveguide into which the nanopillars are embedded. The scale bar in the inset is 10 /xm. (b) DC response: The 
dashed lines indicate the Coulomb staircase. The non-zero slope of the staircase steps is due to thermal broadening. The inset 
gives the circuit diagram with the mechanically tunable tunneling barriers coupling the leads and shuttles. The gate couples 
capacitively with the two electron islands. 



I. RESULTS AND DISCUSSION 



The resulting /^^-characteristic is shown in Fig. 1(b): the application of a DC bias voltage only across the coupled 
pendula leads to an ohmic response, modulated by Coulomb blockade steps. The gate electrode is grounded in these 
measurements. The total Coulomb energy Ec marks the energy required to transport a single charge through an 
island [5T]. From the data in Fig. 1(b) (marked by the arrow) we find a total energy of Eq w 40 meV. The capacitance 
of each island of radius r^p can be estimated to be Cnp — 47reornp(l + (rnp/d)^), where the inter-nanopillar distance is 
d = 80 nm. In accordance with the optimized values that fit the experimental model (see Methods section), we obtain 
Cnpi — 3.2 aF and C'np2 — 3.9 aF, where we used for the islands' radii the values rnpi — 25 nm and rnp2 = 31 nm. 
Since the two islands are placed in series, the total capacitance is w C'~p\ -|- C~p2- Hence, we obtain a total 
capacitance of C ~ 1.76 aF, from which we extract a Coulomb energy of Ec = e^/2C = 41 meV. This is well above 
the corresponding room temperature energy of 26 meV. We stress that room temperature CB for a coupled shuttle 
has not been observed before. Conventional single electron shuttles reveal CB only at low temperature [HI [53], since 
suppression of co-tunneling is not as pronounced as for the coupled shuttles considered in this work. 

In Fig. 2 the response to a pure DC-bias Vdc and gate voltage Vg of the coupled shuttles is studied. /V"-curves at 
altered gate voltages are shown in Fig. 2(a-c). The dashed lines are linear extrapolations of the initial conductance 
slopes and underline the steps caused by single electron charging. Evidently, a typical Coulomb staircase is super- 
imposed on an ohmic response, due to thermal broadening and shuttling [16j. Note that the mechanical motion of 
the shuttles allow exchange of electrons with one electrode while suppressing it on other electrode. Since the applied 
gate voltage changes the potential around the shuttles, the charging energy required for an additional electron on the 
island can be varied accordingly. The trace at Vq = 25 mV in Fig. 2(b) shows an ohmic response without the CB 
steps. As the gate voltage decreases, the Coulomb staircase reappears (Vq = —55 mV), as shown in Fig. 2(c). The 
traces are not symmetric about zero bias, which is due to the slightly different diameters of the nanopillar islands 
and the possible polarization of an oxide layer on the islands [24]. The full DC output signal from a measurement 
run is summarized with respect to both Vq and l^c in Fig. 2(d), revealing Coulomb diamonds. We note that the 
diamond structure compares to that of a single island, due to the common Vg and the similar size of the islands. 
The lack of symmetry of the diamonds with respect to Vdc = is due to the actual difference in island radii, and to 
the difference in capacitance couplings to the source or drain. Although the obtained IV curves do not show clear 
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FIG. 2: Coulomb blockade in the coupled electron shuttle: (a-c) Current Ids traces vs. bias voltage Vdc for gate voltages 55, 
25, and -55 mV. Experimental (d) and theoretical (e) Coulomb diamonds traced in the normalized current, |-^ds|/^ds'', i^i color 
scale representation. The lower borders of the CB regions are represented in red with the Coulomb plateaus depicted in green 
(see color scale). The borders of the CB determined from the theory plots are marked by dashed black lines as a guide to the 
eye, both in (d) and (e). The theoretically derived lines trace the measurements closely. The horizontal lines in (d) indicate 
the previously shown line plots in Fig. lb and Fig. 2(a-c). 



steps due to the thermal environment, the periodic variation of the source-drain offset voltages is seen. The line plots 
for Vg = 55,25, and —55 mV are indicated by the white horizontal lines in Fig. 2(d). To validate the experimental 
results, we compute the current using numerical methods (see the Methods section). The theoretical results shown 
in Fig. 2(e) reproduces the shape and size of the measured Coulomb diamonds and agrees with the estimated size of 
the islands. 

Next we apply an RF signal with zero DC bias in order to tune the mechanical motion of the system, as shown in 
Fig. 3(a). The ability to 'dial in' the shuttling frequency is the fundamentally new feature of this device, as compared 
to classical single electron transistors. We can now trace the DC bias voltages at the three different mechanical modes 
marked in Fig. 3(a), plotted in Fig. 3(b), (c) and (d), respectively. These correspond to commensurate mode numbers 
p/q of the fundamental mode (for more details, see [20]): p/q is 1/6 (circle, / = 87MHz), 5/9 (square, / = 289MHz) 
and 1/1 (triangle, / — 500MHz). The current offset /sb due to spontaneous symmetry breaking (SB) in coupled 
electron shuttles [5D] is subtracted. Dashed lines extend the inclined plateaus to underline the staircase current. 
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FIG. 3: (a) Full frequency sweep of the direct current (Ids) through the coupled shuttle revealing the mechanical mode 
structure for Vdc ~ 0. The colored symbols indicate three distinct mechanical frequencies plotted in (b-d). The lower traces 
indicate the DC signal at those different modes: / = 87, 285 and 500 MHz, respectively. 



Fig. 4(a) exemplifies the dissipative regime: In this limit, a saturation current (frequency independent) is achieved, 
as the effective contact time of the islands with the electrodes is proportional to the inverse of the frequency, and at 
the same time, the current is proportional to the time contact and to the frequency (see Methods section). We observe 
a saturation value for the current of about 2 pA. The measured current steps (black dots) fall on the expected shuttle 
current (red solid line), as predicted from [2] We have to note that their calculation only considered a single shuttle. 
The inset in Fig. 4(a) shows a magnified view of the expected current in the low frequency region. The coupled 
shuttles are operated in the high frequency, dissipative limit (red region). We estimated a;max ~ 5-8 nm, d ^ 17 nm, 
and A ^ 1 nm, giving to ~ 1/(20;). We also have t ~ xlO~^s. Thus, we find the transition frequency to be of the 
order of 5 MHz, which agrees very well with our data. The standard deviation a of the experimental values are found 
to be much lower than the theoretically expected fluctuations, as seen in Fig. 4(b). This can be ascribed to the fact 
that the model assumed a single shuttle unlike the actual device of two coupled charge shuttles in series. However, 
charge fluctuations are even more strongly suppressed off by the mechanical resonances for two coupled shuttles. 

We consider now the conductance level within the plateau. Fig. 4(c) shows two /^-traces at 285 MHz (red squares) 
and 500 MHz (blue circles). The dashed lines give the slopes of the corresponding plateaus in the neighboring 
Coulomb regions. The slope, or differential conductance G — A/ds/AVdc, increases hnearly towards higher shuttling 
frequencies, a consequence of self-excitation. This occurs around Vbc = 0, where a mechanical mode was located, 
and at Vdq > Vc- In the CB regime, the number of transferred electrons per period is quantized in units of e, Ne oc 
(VocC'/e-l- l/2)e. Moreover, in the self-oscillatory regime (above a critical voltage, Vc [H], or around Vdc = 0, where 
the mechanical mode was excited), the charge exchange with the leads takes a fixed quantized value corresponding to 
the thermal equilibrium with the nearest lead, and is exponentially suppressed with the far lead. Following Gorelik 
et al. we estimate the critical voltage Vc to be of the order of the CB voltage, Vc ^ e/C [T5]. In this regime, the 




FIG. 4: Evaluation of time dependent shuttling: (a) The current steps between plateaus A/ds are plotted vs. the mechanical 
frequency / at a fixed gate voltage (Vq = V). The measured current steps (black dots) fall on the expected flat-line, as 
predicted [18]. The inset shows a theoretical plot of the expected current in the low frequency region. The coupled shuttles 
are operated above the cutoff frequency (red region), (b) The standard deviation a of the measured current (black dots). 
This is compared to the theoretical expression for charge fluctuation in a single shuttle (red solid line) as obtained from |18) . 
The strongly reduced co-tunneling leads to a suppression of fluctuations, allowing the observation of CB. (c) JV^-traces for two 
mechanical modes, as seen the slopes in the plateau increase from a shuttling frequency of 285 MHz (red squares) to 500 MHz 
(blue circles), (d) Summary of the differential conductance within the Coulomb-plateau for the different shuttling modes. As 
expected the effective differential conductance is increased to the mechanically enhanced escape rate of the electrons. 



conventional linear-in-frequency relation applies, I ~ 2Nef . Note that for intermediate Fdc values, 0< Vdc < 14, the 
amplitude of the self-oscillations is suppressed. In this region, [2] applies, rather than the linear-in-frequency relation 
(see also the Methods section). The overall results are summarized in Fig. 4(d), where we show the step differential 
conductance in the plateau, A/ds/AVdc, or the slope of the dashed lines of Fig. 4(c). The linearity in frequency of 
the plateau differential conductance G is evident. This demonstrates the self-excitation of the nanopillars, already 
shown to cause ohmic behavior [I3J [2Z1 • 
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II. CONCLUSION 



In summary we observe Coulomb blockade in a coupled electron shuttle operating at room temperature. Mechanical 
clocking of electron transport has been demonstrated by using RF signals. Operating in the high frequency regime 
leads to a strong suppression of current fluctuations. Two regimes of operation are observed in the high frequency 
limit as a function of bias voltage, corroborating theoretical predictions: the dissipative regime, where the current is 
frequency independent, and the self oscillatory CB regime, showing a linear- in-frequency conductance. 



III. METHODS 



To validate the experimental results obtained in Fig. 2(a-d), we compute the current within the orthodox model 
of CB. The accessible states are described by their respective probabilities Pni,n2i where rii is the number of excess 
electrons on the islands i = 1, 2. In the absence of an RF signal, the system can be assumed to be in the dissipative 
regime, where the amplitude of the mechanical oscillations is small and the tunneling current prevails over the 
shuttle current [16]. Changes in the resistance of a particular junction are proportional to e^'^'*', with A being the 
typical tunneling distance, Xi the displacement, and < A in the dissipative regime. In this limit, we suppose l/F 
to be much smaller than the typical time scale of the displacements, and the time evolution can be described in 
terms of coupled master equations !25j , which we solve by exact diagonalization. In the CB limit, we observe that 
d^DsD/^Ds" Sni.na ^«i:n2 ( r^ni,n2 " ^ni.ris): wlicrc (. . .) dcuotcs an eusemblc average. ~f ni,n2 gives the rate of 
tunneling from left to right at one of the junctions for a given configuration {ni,n2\'. 

"^'"^ ^e2i?y exp[AS(ni,n2)/fci3r] - 1' ^' 

with R being the average resistance of the junction and the chemical potential. The theoretical results are shown 
in Fig. 2(e): the CB region below the first transition is shown in red. We stress that a finite conductance is observed 
in this region due to thermal broadening and direct electron shuttling. In the dissipative regime, however, we assume 
temperature broadening to be the dominant source of the ohmic response. 

We now focus on the ohmic response between the differential conductance plateau or current step, A/ds- Con- 
ventionally, one would expect a simple linear relation between this ohmic current and the shuttling frequency of the 
form, / = 2e/ or AJds ~ /, which is typically applied for electrometry purposes |26j . This relation, however, is only 
valid in the limit of large tunneling rates, Fij- ^ with '^ij defined in [l] For our setup, the relation applies 
only in the self-oscillatory regime, achieved in the plateau or around Vdc — 0, where a mechanical mode is excited. 
Our theoretical results indicate that in the inter-plateau region Tij < 5-15 MHz, which is slower than the mechanical 
motion. Hence, we can state that the coupled shuttles operates in the dissipative regime. According to Weiss & 
Zwerger ^8j, we can calculate in this regime the average number of transferred electrons (N) per period: 

^ ' (l + a)(l + l/2a + a2)' 



where a = exp(— io/T), and 2io is the effective contact time, to = u ^ \J^^^{^ + The length Xmax indicates 
the maximum displacement of the mechanical oscillation, and w = 27r/. In the regime of high RF frequencies, the 
effective contact time is shorter than the charge relaxation time, <C t. Thus the number of transferred electrons 
per period is inversely proportional to the frequency of the shuttle oscillation, {N) oc with cx 1//. Then, the 
electron shuttle current achieves a frequency independent saturation value, / = {N)ef cx A/og. 
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